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ABSTRACT 


Theoretical  calculations  of  the  phenomenology  of  an  atmospheric  high-explosive 
detonation  are  presented.  The  change  was  a  100-ton  sphere  (radius,  240.5  cm)  of 
TNT  (loading  density,  1.56  gms/cm3)  whose  center  w=»s  at  an  altitude  of  646.405 
meters.  The  ground,  upon  which  the  spherical  charge  rested,  was  644  meters  above 
se»  it^el.  The  uuu.ei.lcaL  caic»  Lai  ions  lateen  ool  to  ^  "•'■onds  were  pc.r*\;  r"»»d  »»> 
the  CDC  6600  digital  computer  using  SHELL2,  a  two-material  (version  of  the 
SHELL-OIL  code),  two-dimensional  pure  Euierian  hydrodynamic  code.  Air  and  the 
detonation  products  of  TNT  were  the  two  materials  considered  in  the  calculation. 
The  analytic,  self-similar  solution  for  the  detonation  wave  in  TNT  provided  the 
initial  conditions.  Included  are  pressure  and  density  contours,  velocity  vector 
plots,  and  wave  forms  for  19  test  stations.  This  calculation  is  a  representa¬ 
tion  of  the  air  blast  of  Event  6  of  the  DISTANT  PLAIN  test  series  to  be  fired  in 
Canada,  July  1967. 

(Distribution  Limitation  Statement  To.  2) 
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SECTION  I 
INTRODUCTION 

The  United  Spates  is  participating  with  the  United  Kingdom  and  Canada  in  a 
non-nuclear  test  series  taking  place  in  Alherta,  Canada.  The  purpose  of  this 
series,  called  Operation  DISTANT  PLAIN,  is  to  obtain  experimental  data,  using 
high  explosives  and  detonable  gas  o,  which  may  be  applied  toward  solution  of 
problems  concerning  Llart  eff?:":  -.ilifary 

At  the  request  of  the  Defense  Atomic  Support  Agency  (DASA) ,  the  Air  Force 
Weapons  Laboratory  (AFWL)  performed  several  one-  and  two-dimensional  theoretical 
pre-shot  calculations  of  events  held  in  the  summer  of  1966  (Ref.  1).  These  cal¬ 
culations  were  performed  on  the  CDC  6600  electronic  digital  computer  using  SAP, 
a  one-dimensional  hydrodynamic  code  and  SHELL,  a  two-dimensional,  pure  Eulerian 
hydrodynar-’  3  code.  Included  in  these  numerical  experiments  was  the  calculation 
of  such  physical  phenomena  as  the  detonation  process,  that  is,  the  "burn"  of 
the  explosive  substance,  shock  formation  and  growth,  shock  reflection,  mach  stem 
formation,  and  triple  point  path  as  well  as  blast  parameters  such  as  peak  over¬ 
pressures,  dynamic  pressures,  positive  phase  durations,  positive  impulses,  and 
shock  arrival  times  of  the  air  blast  wave.  Agreement  between  theoretical  results 
as  computed  by  the  hydrodynamic  computer  codes  and  the  experimental  data  is 
excellent  (Ref.  2). 

DASA  has  again  requested  the  AFWL  to  perform  theoretical  calculations  of 
DISTANT  PLAIN  Event  6  scheduled  for  July  1967.  This  event  involves  the  detonation 
of  r>  spherical  charge  of  100  tons  of  TNi  resting  on  the  ground.  Reported  in  this 
document  are  the  results  of  these  calculations. 
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SECTION  II 
THE  CALCULATION 

We  performed  the  numerical  calculation  reported  here  on  the  CDC—  .*.CT 
electronic  digital  computer.  The  computer  time  required  to  carry  the  calculation 
out  to  6  seconds  was  approximately  50  hours.  SHELL2  was  the  code  used  to  make 
this  calculation. 

1.  SHELL2 

SHELL2  is  the  two-material  version  of  SHELI -OIL,  a  two  dimensional,  axially 
symmetric,  pure  Eul^rian  hydrodynamic  code.  Since  SHELL-OIL  is  described  in 
considerable  detail  in  reference  1  we  shall  only  present  a  brief  description  of 
Its  computing  method  here.-  However,  we  shall  give  a  more  detailed  description 
of  the  treatment  of  two  materials. 

a.  SHELL-OIL 

SHELL-OIL  solves  problems  In  compressible  fluid  dynamics  by  dividing  the 
region  occupied  by  the  fluid  Into  a  mesh  of  fixed  cells  and  then  forming  finite 
difference  analogs  to  the  Eulerian  hydrodynamic  equations.  The  fluid  la  described 
at  any  Instant  of  time  by  specifying  the  velocity,  density,  and  specific  internal 
energy  of  the  fluid  contained  in  each  cell.  These  values  are  considered  known 
at  the  center  of  each  cell  and  constant  over  It. 

SHELL-OIL  is  a  time-marching  code.  Using  the  finite  difference  analogs, 
it  numerically  integrates  the  hydrodynamic  equations,  updating  the  values  of 
the  thermodynamic  variables  describing  the  fluid  in  each  cell  of  the  computa¬ 
tional  mesh  to  time,  t  +  At  based  upon  their  respective  values  at  time,  t.  During 
each  time  cycle  of  calculation  it  performs  the  following  operations. 

(1)  It  calculates  a  time  step,  At,  such  that 

t  - 1« 

where  the  quotient  (Ax/v)  is  the  minimum  of  the  entire  mesh.  The  velocity,  v, 
may  be  either  a  particle  velocity  or  the  sound  speed  of  a  cell  and  Ax  its  minimum 
dimension. 
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(R)  An  equation  of  state,  p  p(p,J),  assigns  a  pressure  to  each  ceii 

wiicre 

p  5  the  density  of  the  cell  in  gms/cm3 

I  £  the  specific  internal  energy  of  the  cell  in  ergs/gm. 

Vi)  Ph-se  1.  Solution  of  the  equations  for  the  conservation  of  momentum 
and  energy,  neglecting  fluid  motion  (i.e.,  the  convective  terms,  ^  •  V  and 
•  VI,  respectively,  are  dropped)  gives  tentative  new  values  for  velocity  and 
specific  internal  energy  respectively.  That  is.  Phase  1  calculates  only  the  pres¬ 
sure  contribution  to  A^  and  AI.  The  AI  calculation  uses  a  time-centered  velocity, 
1/2  (old  v  +  tentative  new  v)  (Ref.  3). 

(4)  Phase  2.  Solution  of  the  equation  for  conservation  of  mass  gives 
the  mass  flow  across  each  cell  boundary.  Donor  cell  densities  and  velocities, 
interpolated  between  the  donor  and  receiver  cell,  determine  the  amount  of  mass 
flow.  The  velocities  used  are  the  tentative  new  velocities  calculated  in  Phase  I. 

(5)  Repartitioning.  Conservation  of  mass,  momentum,  and  total  energy 
determines  the  final  values  for  mass,  velocity  and  specific  internal  energy  of 
each  cell.  Whenever  some  mass  leaves  one  cell  for  another,  it  carries  with  it 
momentum  (Amu  -  radial  momentum,  Amv  -  axial  momentum) ,  kinetic  energy 

[l/2  Am(u2  +  v2)J  and  internal  energy  (Ami)  where  Am  is  the  amount  of  mass  trans¬ 
ferred  and  u,  v,  and  I  are  the  radial  velocity,  axial  velocity,  and  specific 
internal  energy  of  the  donor  cell.  These  values  are  subtracted  from  the  donor 
cell  and  added  to  the  receiver  cell.  Conserving  mass  gives  the  final  value  of 
mass  for  each  cell.  Conserving  momentum  gives  the  final  values  of  velocity. 
Conserving  total  energy  gives  the  final  value  of  specific  internal  energy. 

b.  Two-Material  Considerations 

The  first  requirement  necessary  for  taking  a  second  material  into  consi¬ 
deration  in  the  hydrodynamic  calculations  is  additional  computer  storage.  Storage 
is  required  for  a  second  mass  block  and  a  second  internal  energy  block;  that  is, 
the  amount  of  mass  and  internal  energy  of  the  second  material  present  in  each 
cell  must  be  known.  The  operations  of  SHELL2  proceed  as  follows. 

(1)  The  calculation  of  a  time  step  in  the  two-material  version  of  the 
code  is  exactly  the  same  as  for  the  one-material  version  outlined  above. 

(2)  The  appropriate  equation  of  state  assigns  a  pressure  to  each  cell 
containing  only  one  kind  of  materia).  For  those  cells  which  contain  a  mixture 
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of  both  materlalr,  the  requirement  of  pressure  continuity  across  a  materia] 
interface  defensives  the  procure  o.  thecn  cells.  That  is 


P  *  Pi  (*>,  >li)  m  \  ^2»T2) 


where  Ij  and  I2  art  the  oecific  internal  energies  of  the  first  and  second 
material,  respectively,  ir  a  nixed  cell.  The  densities  p'j  and  p'  are  the 
"effective"  densities  of  each  material: 


P 


1 

2 


Ml 

V' 

V1 


v  *  v'+  V* 

1  2 

where  Mj  and  M2  are  the  masses  of  each  material  contained  ir  a  mixed  cell  whose 
volume  is  V.  Since  all  quantities  are  known  except  the  "effective"  volumes  of 
each  material,  the  pressure  of  a  mixed  cell  must  be  determined  by  iteration. 

(3]  Fn.'  .  1  of  the  two-material  calculation  proceeds  exactly  in  the  same 
manner  as  in  the  one-material  case  outlined  above.  The  only  exception  is  that 
for  cells  containing  a  mixture  of  both  materials,  the  change  in  specific  internal 
energy,  AI,  must  be  distributed  between  the  two  materials.  Assuming  that  changes 
ir.  a  mixed  cell  take  piece  adiabatically ,  the  change  in  internal  onergy  of  each 
material,  AI^  and  AI2,  respectively,  is 


AI, 


VA 

— —  M,  I. 

V  1  V 


AI 


*1-1 


Mi  Ti  12 


Y^l 

r2 


Y„-l 

— M,  I, 

y,  2  2 


AI 


AI 


Y,-l 


y2-i 


„i  +_  m2  I 


Yl  1  1  Y2 


(T^f ,  3)  where  the  y's  are  the  ratios  of  specific  heats  of  each  mater -al, 
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<A)  Phase  2  also  proceeds  in  exactly  the  same  manner  as  in  the  one- 
material  case.  Here  again  the  exception  is  that  the  mass  transported  from  a 
mixed  cell  must  he  divided  between  the  too  materials.  The  apportionment  is  pro¬ 
portional  to  the  amcr.it  of  mass  of  each  material  present  in  the  donor  cell. 


Since  Eulerian  cells  are  homogeneous,  which  leads  to  a  false  material 
diffusion,  special  care  must  be  hafcen  not  to  allow  one  material  to  erroneously 
diffuse  through  the  otner.  To  counteract  this  possibility,  before  a  given 
material  is  allowed  to  transfer  to  a  cell  which  does  not  have  any  material  of 
this  type,  its  presence  in  the  donor  cell  must  be  in  sufficient  quantity  such  that 
it  is  at  least  a  predetermined  minimum  percentage  of  the  total  mass  of  the  donor 
cell. 


(5)  Repartitioning  i r  also  similar  tc  that  in  the  one-material  version. 
For  cells  containing  one  material  it  is  exactly  the  same.  For  mixed  cells  con¬ 
servation  of  Momentum  is  also  the  same: 


u 


(n+1) 


p(n-Kl)  _ 

(n+J)  (n+1) 

1  +  ll2 


_(n+l) 


z(n+l) 


where  is  the  pest-transport  radial  momentum  in  a  mixed  cell  and  Z 

is  the  post-transport  axial  momentum  in  a  mixed  cell.  The  updated  values  of 
specific  internal  energy  for  each  material  in  a  mixed  zone  come  from  the  inde¬ 
pendent  conservation  of  total  energy  for  each  material. 
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(n+l)  _  TE2 


(n+l) 


(n+l) 


M, 


-  1/2  [(un+1)2  +  (vn+i)2] 


where  the  TEs  are  the  post-transport  total  energies  of  each  material  In  a  mixh^ 
rone. 

2.  Initial  Conditions 


The  SHELL2  computational  mesh  contained  123  cells  in  the  z-direction  and  121 
cells  in  the  radial  direction  for  a  total  of  14,883  cells.  All  cells  except  the 
bottom  row  of  121  cells  comprised  the  active  mesh.  Each  cell  in  the  bottom  row 
was  544  meters  high  (Az)  and  5  cm  wide  (Ar) .  Each  cell  in  the  active  mesh  was  a 
5-cm  square.  Therefore,  the  overall  nesh  extended  from  sea  level  to  an  altitude 
of  650.1  meters  and  out  to  a  radius  of  6.05  meters;  and  the  active  portion  of  the 
mesh,  from  644  meters  (ground  level)  to  an  altitude  of  650.1  meters.  The  bottom 
boundary  (ground  level  at  644  meters)  was  reflective;  the  top  and  right  boundaries 
were  transmittiva  and  thus  subject  to  rezone  (Ref.  1)  in  their  respective  direc¬ 
tions.  The  left  boundary  (r=0)  was  the  axis  of  symmetry. 


All  calls  in  the  active  mesh  with  the  exception  of  those  contained  in  a 

sphere  whose  center  was  on  the  z-axis  at  an  altitude  of  646.405  meters  and  whose 

* 

radius  was  240.5  cm  contained  an  ambient  (real)  atmosphere  (Ref,  1).  A  similarity 
solution  of  the  flow  field  behind  the  detonation  front  provided  values  of  velocity, 
density,  and  specific  internal  energy  which  were  fitted  into  those  cells  contained 
in  the  sphere  (see  figure  1),  Table  I  contains  the  data  of  the  detonation  wave. 


A  Lagrangian  interface  of  715  trace  particles  simulated  the  motion  of  the 
TNT-air  interface.  These  particles  were  initially  placed  in  -  'semicircle  whose 
center  was  on  the  z-axls  at  altitude  of  646.405  meters  and  whose  radius  was 
240.5  cm.  There  were  also  701  more  trace  particles  placed  in  those  cells  con¬ 
taining  the  detonation  products  of  TNT.  These  particles  did  not  affect  the  hydro- 
dynamic  calculation  in  any  way.  Their  motion,  calculated  with  interpolated  cell 
velocities,  depending  on  their  position  in  a  cell,  merely  indicated  the  material 
interface  and  gave  a  general  picture  of  the  motion  of  the  TNT  detonation  products. 

Also  introduced  into  the  calculation  were  19  "test  stations"  where  the 
overpressures,  dynamic  pressures,  impulses,  and  velocities  were  calculated  each 
cycle.  The  positions  of  these  test  stations  were  chosen  to  coincide  with  the 
position  of  the  experimental  stations.  See  table  II  and  figure  2. 
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The  two  materials  considered  in  this  calculation  were  the  TNT  detonation 
products  and  air.  The  corresponding  equations  of  state  used  in  the  calculation 
were  the  LS^.K  equation  of  state  for  the  former  and  the  Doan-Nickel  equation  of 
state  for  the  latter. 

Table  I 

DETONATION  WAVE  FOR  TNT  (p  -  1.56  grns/cc) 

Specific 

r  Internal  Density  Velocity 

(cm)  Energy  (ergs/gm)  (gras/cc)  (cm/ sec) 

0.  2*703638+10  1 «  2374  CE  +  00  0. 

1.09026E+0?  2.703638  +10  1.23740F+00  0. 

1.09788E+0?  2.70410E+1O  1.23756E+00  3*96^308  +  01 

1 • 1 5627E+0?  2*  73099E+ 1 0  1  .24674E  +  00  2*289048  +  03 

1.22050E+0?  2.77351E+10  1*261168+00  5.697138+03 

1.27617F+0?  2.01772P+1O  1 *276048+00  9. i 0*2?- +0 ^ 

1.32733E+0?  2.96339F+10  1.29128E+00  1.251338+04 

1 • 3755 1 E+02  2*91 039F+ 1 0  1 .30685E+00  1.59214E+C4 

U42146E+0?  2.95862E+10  1.32269E  +  00  l  . 932'->5F  +  04 

1.46564E+0?  3* 00802E+ 1 0  1.33879E+00  2*273768  +  04 

1.S0832E+0?  3*0585 1 E+ 1 0  1.35510E+00  2. 614578+04 

1.54968E+0?  3*110058+10  1.37162E+C0  2. 985338+04 

1.5J966E+0?  3*1 6260E+ 1 0  1.38833E+00  3. 296 1 M  +C 4 

1. 628948+0?  3.2161 JE+1C  1.40520E+00  3.637Cr-_  +  04 

l • 66698E  +  0 2  3.27055E+10  '.42222E+00  3.9'7918+04 

1.704C2E+C?  3.32587E+1C  lt*3937E+00  4 , 3  -c* 

1.7400LE+0?  3*382048+10  1 *45665€*00  4.659438*04 

1  •  775  1  9E  +  0?  3.43OQ3P+10  1.47403E+0C  ,  OOC.'*  .  ■‘•''6 

1 *e093«E+02  3.49680E+1C  1.4U?*j8+oo  5.341048+04 

l*S42*tet+02  3.55532F* 1 0  I *50«07E+CC  5.68»8f?*C4 

U874*4€+0?  3.614558+10  I  .526708+00  6 . 022668 * 0 4 
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T  (eonc'd) 


Table 


Specif ic 

K  Internal  Density  Velocity 

(cm)  Energy  (ergs/gm)  (gms/cc)  (cm/sec) 


j  ,9061  6fc'+nr>  3 ,  r,',  4  4  5F  +  I  n 
1.936S3E+0?  3".  7350  IF*  10 

1  .  96594E  +  0/#’  3.  79%  1 7F+1 0 

1SQ94  39E40(2  3.8%702f- +  i  0 
2.02186E+0?  3.92020F+10 

2  «  04835F  +  02  3.99200F+tO 
2  •  0  7  385E  +  02..  4.04626E+10 
2.09837F+02  4.  109O6E*  1  0 
2.121 88F+O9  4 • 1 7407F+ 1 0 
2.14439E  +  0^>  4.23854F+1C 
2.165O0E+0?  4.30334E+10 
2*1 864 1 E+O?  4.36044E+10 
2.20592E+02  4.43380E+10 
2.22443E+0?  4.49939E+10 
2.24195E+0?  4.56517F+10 
2.25848E  +  02  4.631  lPF+m 
2.27403E+0?  4  *6971 6E+ 1 0 
2.28862E+0?  4.76330^+10 
2.30225E+G3  4.P2Q50F+I0 
2.314O3E+0?  4.qQR7lF+10 
2  •  3266°E-**0  ?  4.o61r,2F+!0 
2 • 33764E  +  0')  8.02M3flF+ l 0 


1  .S443OE+0C  6.3634 *7F  +  o  + 

1.56213E+00  6.7042PF+04 
1  *F700  1  E  +  00  7.045008+04 
I.59771E+00  7.385«OF+04 
1.61 553E+00  7.72671E+04 
1  .63336F+00  8 . 06 7026 +0  4 
1.651 18E+00  8.40833E+04 
1  .66900E+00  8.74O14E+04 
1  .6R679F+C0  9.0B995F+04 
1.70454E+00  9.43076E+04 
1 .72226E+00  9.771F7E+04 
1 .73994E+00  1.011248+05 
1 .75755E+00  1 .04532E+05 
1.77510E+00  1 .070408+05 
1 .79258E+00  1 . 1 1 348E+CF 
1  .8noQpF-j.no  1  ■  1 475f»F 4  n« 

1 .82729E+00  I.1S164E+0F 
1 .84451F+00  1.21F7?r+0r 
1.86163E+00  1.24O80r+nr 
1  .87B65E+00  1»2S3E«Z40* 

1 .8«555E+00  i.317«7E+0B 
1  .ni  233E  +  C0  1  .3spr?E+0*v 
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Table  I  (cont'd) 


Specific 

R  Internal  Density  Velocity 

(cm)  Energy  (ergs/gm)  (gma/cc)  (cm/sec) 


2.34749F+&2  «.09416F+10 
2.356F6EHD?  5.16014E+1C 
2  a 36478E+0?  F,225«PF+t 0 
2 • 372 ! 6E+02  5.29165F+10 
2.37872E+0?  5.3571 3E+10 
2.38450E+0?  S.42238E+10 
2.38950E+02  S.48738E+10 
2*39376E+0?  5*  352  HE*  10 
2*39729E-*-Q:>  5.616S3E+10 

2.40013E+02  3.68063E-M0 
2.40230E+02  3.74438F+10 

2.40381E+0?  5.80775E-MO 
2 #4047 1*6402  S.87072E+1C 
2 • 4QSOOE+ 02  3.93328E-M0 


i  .9?n«eE+no 

1  .  7=161  ',tr  +  C^ 

1  .94551E+00 

1 .4202 1 F +05 

1  .96191E+00 

1  ,45429E+0~ 

l .97B16E+00 

1 . 4QO37E  +  0c: 

l .99427E+00 

1 *S2245E+0F 

2.01023E+00 

1 • 55653F  +  05 

2.02604E+00 

1 .59061E+05 

2.04170E4-00 

1 .6246OE+03 

2.05720E+00 

1  t6587P.F+03 

2.07253E+00 

1 .69266E+05 

2.08rtt}£4C!0 

*  • . 

1 .726948 +05 

2.1«27t»£+00 

1 • 76 1 02E+05 

2*1 1793(400 

1.79S10E+05 

oo 

1 .829 1 8E+05 
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Table  II 

STATION  POSITIONS 


Station 

number 

Range 

(ft) 

Altitude 

(ft) 

1 

25 

0 

2 

35 

0 

3 

48 

0 

A 

51 

0 

5 

69.3 

0 

6 

85 

0 

7 

115 

0 

8 

170 

0 

9 

232.4 

0 

10 

270 

0 

11 

340 

0 

12 

400 

0 

13 

550 

0 

14 

850 

0 

15 

0 

25 

16 

0 

35 

17 

12 

0 

18 

15 

0 

19 

21 

0 

<»)  aoruuiv  (u)  saruinv 


1? 


Figure  2.  Position  of  Test  Stations. 
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SECTION  III 
RESULTS 

The  SHELL2  calculation  began  with  an  analytic  (Taylor)  solution  for  the 
detonation  as  the  starting  conditions.  The  detonat  on  wave,  which  reached  the 
surface  of  the  TNT  (p^  »  1.56  grams/cm3)  sphere  at  0.35285  millisecond  following 
initiation  of  detonation,  had  the  following  parameters: 

Detonation  wave  velocity  6,8.16  *  105  cm/sec 

Pressure  1.945  *  101 1  dynes/cm 

Density  2.132  gms/cm3 

Particle  velocity  1.829  *  105  cm/sec 

Specific  internal  energy  5.933  *  1Q1C  ergs/gm 

Figures  3  through  9  give  a  summary  of  blast  narameters  obtained  from  this 
calculation.  The  appendixes  contain  plots  giving  a  detailed  history  of  the 
calculation. 

1.  grec-FlelJ  riienomenclogy 

The  calcination  began  with  the  free  expansion  of  the  detonation  products  into 
air.  The  ground,  which  SHELL2  treated  as  a  perfectly  reflecting  plane,  destroyed 
the  symmetry  of  the  problem  at  the  very  beginning.  The  effects  of  the  ground 
were  first  evident  at  the  point  of  the  tangency  where  the  TNT  sphere  was  resting 
on  the  ground.  As  the  calculation  progressed,  the  ground  effects  spread  very 
quickly  over  the  entire  region  of  interest. 

The  phenomenology  In  the  vertical  direction  from  the  point  of  burst  resembled 
that  of  free  air  longer  than  in  any  other  direction.  When  the  detonation  wave 
reached  the  s '.”’race  of  the  TNT  sphere,  the  pressure  profile  had  a  peak  of 
1.943  x  10'*  dynes/car  at  the  surface  <r*Z40.5  cm)  and  decreased  to  4.47  •  10- v 
dynes/car  at  r  *  110  cm.  This  pressure  remained  constant  back  to  r  *  0.  At 
this  time  the  detonation  products  escaped  into  the  air  with  maximum  escape  velo¬ 
cities  of  5  «  10^  cta/sec.  This  rjpid  expansion  of  the  detonation  products 
accelerated  tit*  air  Just  in  front  of  it,  thus  producing  a  shock  wave.  Meanwhile 
a  second  wave,  a  rarefaction  wave,  propagated  into  the  detonation  products, 
decreasing  the  pressure  as  it  proceeded  inward.  See  figure  3.  The  rarefacti-n 
wave  reached  the  center  at  about  0.95  millisecond.  A  second  shock  wave  formed 
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Fi$«r«  3.  Free  Expansion  of  the  Dttocution  ?;o*uc 
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Figure  7.  Dynamic  Pressure  Impulse  vs.  Ground  Range. 
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in  the  detonation  products  Just  inside  the  TNT-air  laterfacu.  It  was  initially 
very  weak  but  gradually  gained  in  strength.  Because  of  the  high  expansion 
velocities,  the  second  shock,  though  following  th.  rarefaction  wave  into  the 
detonation  products,  was  initially  moving  outward  in  space.  As  the  energy  of 
the  detonation  products  dissipated,  the  velocity  of  the  second  shock  became 
larger  than  the  velocity  of  the  detonation  products,  and  so  the  second  shock 
began  tc  move  Inward  toward  the  center.  Thi*  occurred  at  approximately  25  milli¬ 
seconds.  After  this  time  the  phenomenology  no  longer  resembled  that  of  free  air 
as  seen  in  the  plots  of  appendix  I. 

The  phenomenology  for  free  air  is  fully  documented  in  references  1  and  4. 

2.  Effecte  of  the  Ground 

At  the  point  of  tangency  the  TNT  detonation  products  could  not  expand  nor 
could  they  reflect  because  the  mass  above  was  still  moving  outward.  This  caused 
a  build-up  in  density  and  a  correjponding  rise  in  pressure  along  the  ground  with 
the  only  means  of  relief  being  u  pith  outward,  parallel  to  the  surface  o i  the 
ground.  Thi6  motion  and  the  resulting  large  horizontal  velocities  can  be  seen 
in  the  plots  of  append!.'.  I  as  early  as  0.5  trlllisecond.  Along  this  path  the 
detonation  products  expanded  more  rapidly  than  in  any  other  direction,  thus 
creating  a  "toe"  to  an  otherwise  spherically  expanding  shock  wave.  Inis  toe 
became  the  dominating  feature  of  the  calculation  and  remained  prominent  to  a  time 
or  300  milliseconds.  At  this  time  the  shock  had  reached  a  distance  of  approxi¬ 
mately  700  feet.  From  this  time  on,  the  main  shock  front  resembled  that  of  a 
hemispherical  surface  detonation.  By  500  milliseconds  the  toe  was  gone  -md  the 
main  shock  front  was  approaching  a  spherically  symmetric  expansion. 

Meanwhile  the  detonation  products,  driven  by  the  large  pressure-density 
gradients  along  the  ground,  accelerated  outward  until  they  reached  the  flow  field 
immediately  behind  the  out ward -moving,  toe-shaped  shock  (8  millireconds) .  In 
this  region  the  pressure  gradients  of  the  toe  imparted  s  positive  vertical 
velocity  component  to  an  otherwise  horizontal  flow  field.  Hence,  the  detonation 
products,  after  reaching  the  shock,  passed  vertically  alot.g  its  back  edge  into  a 
region  of  smaller  horizontal  velocities.  Other  detonation  products  traveling 
along  the  ground  with  larger  horizontal  velocities  passed  underneath  those 
products  which  were  picked  up.  Therefore,  the  TNT-air  Interface  formed  a  hook 
to  the  rear  of  the  expanding  shock  along  the  ground.  The  structure  of  the  toe 
was  such  that  the  peak  pressure  hsd  its  greatest  radial  extent  along  the  ground, 
and  the  radial  position  of  the  peak  decreased  with  altitude.  Hence  behind  the 
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shock  along  Che  ground  Che  flow  field  had  negaCive  verdcal  velocity  cotnponenCs. 
This  caused  Che  deconadon  produces  which  had  been  picked  up  Co  move  back  down 
again,  thus  forming  a  loop.  The  detonation  products  reached  a  radius  of  approxi¬ 
mately  475  feet  and  a  height  of  approximately  30  feet  at  a  time  of  500  milli¬ 
seconds.  The  entire  motion  of  the  products  can  be  seen  n  the  plots  of  appendix  I. 

3.  The  Second  Shock 

Section  1  above  describes  the  formation  and  movement  of  the  second  shock  in 
the  vertical  direction.  If  this  were  a  detonation  in  completely  free  air,  this 
second  shock  would  converge  on  the  burst  point  and  reflect  there.  It  would 
subsequent!)  travel  outward  to  the  TNT-air  interface  through  which  the  second 
shock  would  be  transmitted  and  a  third  shock  would  form  there  and  travel  inward 
repeating  the  process  of  the  second  shock.  See  reference  4.  In  this  fashion 
a  whole  series  of  shocks  would  be  propagated. 

However,  the  presence  of  the  ground  in  this  calculation  vastly  modified  this 
sort  of  phenomenology.  The  second  shock  front  was  very  asymmetrical  and  became 
more  so  as  time  went  on.  Its  motion  was  also  asymmetrical  and  quite  complicated. 
The  plots  of  appendix  I  can  describe  the  motion  of  the  second  shock  front  best. 

A  simplified  version  is  as  follows:  The  second  shock  starts  its  inward  motion 
at  the  top  near  the  outer  edge  of  the  detonation  products.  This  inward  motion 
of  the  second  shock  spreads  to  the  adjoining  portions  of  the  shock  front.  The 
convergence  of  the  shock  on  the  center  of  burst  is  modified  in  such  a  way  as  to 
cause  the  shock  to  bee  ue  nearly  planar,  pass  vertically  through  the  center  of 
burst,  and  reflect  from  the  ground. 

4.  Effect  of  the  "Toe"  on  Blast  Parameters 

The  blast  parameters  near  the  ground  were  drastically  modified  by  the 
reflecting  surface  and  subsequent  formation  of  the  toe.  These  effects  were  not 
simply  close-in  effects  but  remained  important  throughout  the  entire  region  of 
interest.  Increased  overpressures  were  noted,  but  of  greater  importance  were  the 
much  higher  dynamic  pressures  observed  than  would  be  expected  from  an  equivalent 
yield  hemisphere  (Fig.  10).  The  increased  dynamic  pressures  were  further  compli¬ 
cated  by  large  vertical  components  of  velocity  and  dynamic  pressure,  caused  by 
the  increased  pressure  of  the  material  near  the  ground  accelerating  the  region 
above  it.  This  Is  contrary  to  expectations  for  a  hemispheric^  configuration 
of  tha  shock  which  would  have  negligible  vertical  velocities.  The  vertical 
velocities  near  tha  ground  at  the  shock  front  are  positive  giving  the  jnock  an 
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Inclination  of  over  15  degrees  to  the  vertical,  to  a  radius  of  over  350  feet. 
Somewhat  behind  the  shock  front  Is  a  region  of  negative  vertical  velocities 
causing  a  sudden  reversal  of  the  vertical  flow  at  any  point  within  several  meters 
of  the  surface.  This  effect  can  ba  seen  in  the  vertical  velocity  and  vertical 
dynamic  pressure  station  plots  of  Appendix  II. 

The  toe  did  no  allow  normal  shock  formation  until  It  had  reached  a  radius 
of  approximately  150  feet  and  by  200  feet  the  shock  was  well  formed.  At  about 
150  feet  a  sudden  change  occurred  In  positive  phase  duration  and  was  the  most 
dramatic  indication  of  normal  shock  formation.  At  this  same  radius  similar 
deviations  from  a  smooth  curve  are  noted  in  all  other  variables. 

Beyond  the  200-foot  radius  the  shock  retains  some  of  the  characteristics 
imparted  to  it  by  the  toe  of  the  detonation  products. 

The  parameters  appearing  in  this  report  were  taken  directly  from  the  computer 
output.  No  extrapolation  was  used  in  either  overpressure,  dynamic  pressure,  or 
radius.  Extrapolation  was  considered  neither  important  nor  significant  over  the 
greatest  portion  of  the  region  of  interest.  Extrapolation  procedures  have  been 
used  on  previous  occasions  (Ref.  1)  and  the  effect  of  extrapolation  can  be 
accurately  predicted. 

The  overpressures  to  a  radius  of  300  feet  should  be  within  experimental 
error.  Beyond  this  radius  the  peak  overpressure  drops  somewhat  below  the 
expected  overpressures,  and  by  600  feet  are  about  10  percent  low. 

The  one-dimensional  SAP  calculation  has  sufficient  resolution  that  extrapo¬ 
lation  is  not  necess«j.y.  Figure  10  includes  a  comparison  of  the  SHELL,  vertically 
upward,  and  the  SAP  calculation.  To  a  radius  of  400  feet  no  significant  devia¬ 
tion  is  noced. 
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SECTION  IV 
CONCLUSIONS 

This  calculation  ahovs  that  Asymmetries  in  the  flow  field  as  characterized 
by  the  "toe"  moving  along  the  giound  will  be  a  dominant  feature  of  DISTANT 
PLAIN  Event  6  phenomenology.  These  aaymeaetries  will  be  present  throughout  the 
period  of  interest  (0.5  ms  to  300  ms),  a  much  lunger  time  than  previously 
anticipated  (Ref.  5). 

One  should  also  expect  higher  overpressures,  higher  velocities,  and  higher 
dynamic  pressures  farther  out  along  the  ground  because  of  the  longer  duration  of 
these  asymmetries.  Large  vertical  velocity  and  dynamic  pressure  components 
will  be  associated  with  this  toe  just  above  the  ground  and  extending  in  altitude 
to  the  unperturbed  portion  of  the  shock  front. 

The  calculation  reported  here  treated  the  ground  as  a  perfectly  reflecting 
plane.  However,  this  will  not  entirely  be  the  case  in  the  experiment.  Some 
of  the  blast  energy  will  be  transformed  into  ground  motion  and  crater  formation. 
This  energy  may  be  as  high  as  10  to  15  percent.  So  the  actual  effects  of  the 
ground  will  modify  the  pbenomenolcgy  reported  here.  However,  these  calculations 
provide  an  upper  limit  of  what  can  be  expected. 
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A^ENDIX  I 

PRtdSUKE  AND  DENSITY  CONTOURS  AND  VELOCITY  VEC'CL  PLOTS 

Appendix  I  contains  pressure  and  density  contours  and  velocity  vector  plots 
giving  a  detailed  history  of  the  calculation. 

On  all  graphs,  the  integers  along  the  top  are  the  indices'  of  the  cnoputa- 
tiona]  mesh  in  the  radial  direction,  The-  incegirs  along  the  right  are  the 
indices  in  the  axial  direction.  The  computational  mes  of  123  *  illj. 

(axial  *  radial)  colls  indicates  the  active  portion  of  *“he  me<;h.  it”  altitude 
scale  gives  the  height  above  ground. 

On  all  graphs,  the  dots  represent  the  location  of  the  TNT  detonation 
products.  The  heavily  dotted  curve,  which  in  these  graphs  appears  ns  a  thick, 
solid  curve,  represents  the  TNT-«ir  interface. 

On  tve  pressure  contours,  each  line  connect?  oints  of  constant  pressure. 
Each  line  has  a  number  which  corresponds  to  a  value  of  pressure  in  dynes 'ca2 
given  in  the  contour  scale. 

On  the  density  contour,  each  line  cennec  points  of  constant  density. 

Each  "  ine  has  a  number  rhich  corresponds  to  a  \alue  of  density  in  g*s/ca3  giver, 
i  the  cc  “'hir  scale. 

The  velocity  vectors  in  the  velocity  vector  plots  give  the  magnitude  and 
dir-ci'on  of  the  velocities  present  in  the  field.  The  velocity  scale  appearing 
in.  the  upper  right-hand  corner  of  each,  plot  corresponds  to  the  magnitude  of  each 
vector  in  a  plot.  The  tail  oi  each  vector  i;  located  at  the  center  of  the  cel.. 
The  direction  of  the  am  indicates  the  direction  cf  the  velocity. 
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APPENDIX  II 


WAVE  FORMS  AT  THE  TEST  STATIONS 

Appendix  II  contains  the  tracings  recorded  at  each  test  station  of 
overpreacure,  dynamic  pressure,  overpressure  impulse,  dynamic  pressure  impulse, 
and  velocity.  These  plots  are  made  on  two  different  time  scales,  the  first  being 
terminated  at  the  end  of  the  positive  chase  and  the  second  running  from  just 
before  arrival  time  to  600  ms. 
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1j.  AS3TRACT 

Theoretical  calculations  of  the  phenomenology  of  an  atmospheric  high-exploc ive 
deton;  'ion  ore  presented,  rl.a  charge  was  a  100-ton  sphere  (radius,  2  ♦0.5  cm^  of 
TNT  (loading  density,  1.56  gms/cm3)  whose  center  was  at  an  altituae  of  646.405 
meters.  The  ground,  upon  which  the  spherical  charge  rested,  was  644  meters  above 
level.  The  numerical  calculations  taken  out  to  6  seconds  were  performed  on 
the  CD(.  6600  digital  computer  using  SHELL2 ,  a  two-material  (version  of  the 
SHELL-OIL  code),  two-dimensional  pure  Eule.rian  hydrodynamic  code.  Mr  and  the 
detonation  products  of  TNT  were  the  two  material?  consider -d  io  the  calculation. 
The  analytic,  self-similar  solution  for  the  detonation  wave  in  TNT  provided  the 
initial  conditions.  Included  are  pressure  and  density  contours,  velocity  vector 
plots,  and  wave  forms  for  19  test  stations.  This  calculation  is  a  represen¬ 
tation  of  the  air  blast  of  Event  6  of  the  DISTANT  PLAIN  test  series  tc  be  fired 
in  Canada,  July  1967, 
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